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Abstract 

We present a general study of the dynamical properties of Anisotropic 
Bianchi Universes in the context of Einstein General Relativity. Integra- 
bility results using Kovalevskaya exponents are reported and connected to 
general knowledge about Bianchi dynamics. Finally, dynamics toward sin- 
gularity in Bianchi type vill and ix universes are showed to be equivalent 
in some precise sence. 

1 Homogeneous Universe and Bianchi models 

Considering the usual synchronous frame of General Relativity ^ ds^ — g^,j dx^ dx^ 
Qij dx^ dx^ — dt^ . A Universe is said homogeneous when there exists an isom- 
etry group wich preserves the infinitesimal spacial lenght c?P = cjij dx^ dx^ A 
characterization of the isometry group is possible writing structure constants 

C^,^{d.e'^^-d,e'r)eiel (1) 

where dx' = dy^ . Constants C^^^ are tensorial, low components are antisy- 
metrics, and follow the Jacobi rule : 

y-^ e d , & d , e d r» 

Decomposing C^^, = Eabd N'^'^ + Aa — 5^ At , where iV°^ is a symetric tensor, 
one can show that equivalence classes of homogeneous universes are equivalence 
classes of N"'' with N°'''Ab = 0. Without less of generality, the symetry of 7V°^ 
allows us to write 

jyab ^ 



and Ab = [a, 0, 0] (3) 



ni 

712 

na 

^without contrary indications greek indexes run from to 3, latin indexes from 1 to 3, 
metric signature is (—,+,+,+), e and <5 are respectively the completely antisymetric Levi- 
Cevita tensor and the Kronecker symbol 
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Models split then into Class A with a — and Class B with a ^ and can be 
arranged in the well known Bianchi models 
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2 Einstein Equations 

Following P, one writes ds^ = 7 (r) u'oj^-N"^ (r) dr^ with 7 (r) = diag [e^i(^\ e^^^^^, e'^^'^^)] 
and dt = N (r) dr. The so called invariant differential forms basis are linear 
combinations of dxi with exponential or trigonometric coefhcient in Xi (See PP). 
Finally, r and N (r) are respectivelly the conformal time and the lapse function. 



2.1 BKL Formalism 

This formalism was introduced in the 70's by 3 . Filling Universe by a barotropic 
fluid with pressure P and energy density p such that P = (T — 1) p, taking 
N"^ (r) = V'^ = e^^^^^^^^ for the lapse function, some algebra then gives from 
Einstein equations, the equations for the dynamics of Bianchi Universes 
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Em — ~4xp and a / stands for d/dr 
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2.2 Hamiltonian formalism 

This formalism was introduced by 0], at almost the same time than the prece- 
dent. It consists to diagonalize the quadratic form Ec by introducing new vari- 
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Figure 1: Projection of vectors on (61,62) plane 



ables such that 
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q := [qi 92 93]^ - M [Ai A2 As]^ 



p:=[piP2P3]' =M [A'.A'^A'^Y =q' 
Dynamical equations for Bianchi Universes then becomes 

dH , dH ^ , dH , dH 

: Pi,2 = --^ and 93 = — , p =- — 

dpi, 2 091,2 op3 dqs 



9l,2 



(5) 



(6) 



where H — ^ {p,p) + ^^fc^e with the following products 



Vx,yeM^ (x,y) := +xi?;i+a;2y2 + a;3y3 and (x, y) := -xij/i - 2:22/2 +3:3^3 

(7) 

and the constants ki :~ 2nin2, k2 := 2riin3, :— 2n2r7,3, k^ :— —n\, k^ :— —n\, 
fee := — n§, k-j — —ApoX- The set of vectors a^^i,... ,7 is highly symetric (see Fig. 
^ : It allows to use algebraic techniques based on Lie algebra in some dynamical 
treatments of this problem. 



3 Integrability of Bianchi models 

3.1 Kovalevskaya stuff and autosimilar ODE 

If an ODE 

X = f (x) with X G M" and ' = d/dt 
admits a self-similar solution 
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with the weight vector g and the constant vector c lying respectively Z" and 
R". One can then prove that the linearized system around x, admits too an 
auto-similar solution 



di [t - to) 



, dji (t to) 



'On 



(10) 



where d is a constant vector, and, this time the Kovalevskaya vector k, wich 
components are the so-called Kovalevskaya exponents, lies C". In practice, one 
can compute this vector : it is easy to show that Kovalevskaya exponents are 
eigenvalues oi K := D [f (x)] (c) -I- diag (g). A theorem by Poincare shows that 
each component of the non linear general solution of equation © is on the form 



X^ (t) <x{t~ to) ^' S {t- to) 



.{t-tof 



(11) 



where S [.] stands for a multiple series. This result is at the basis of a sufficient 
condition (known as Yoshida's Theorem see [S],|H]) : if all the Kovalevskaya 
exponents are in Q, then the system is algebraicly integrable. 



3.2 Kowalevskaya exponents and Bianchi Universes 

Such a work was Pioneering by [7j , applied for the first time to some special cases 
by ID, and more recently, developed with some imprecisions and incompletude 
by 0. We present here results obtained in a precise way in [Tn|. 
Introducing new variables 

{u e M^, Ui=i,...,7 := (ai,p) 
(12) 
V e W , t;i=i,...,7 := exp (ai,q) 

Einstein's Hamiltonian equations become 

V[ = UiVi 

Vf = 1, • • • , 7 <( 7 with W^J := -kj (a,, a^) (13) 

This dynamical system admits an autosimilar solution x = [At^^, /ii^'^] where 
the constants [A, /i] e x are solutions of the algebraic system 

Wij fij = -Ai 

V* = l,---,7 { (14) 

Ai = — 

Taking into account that Rank(VF) = 3, it exists 45 distincts non trivial 
solution for system H14() , and then 45 sets of 14 Kovalevskaya exponents for all 
Bianchi Universes. Analysis of such sets let us claim that : 
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• In vacuum or with stiff matter (F = 2), excepted Bjx and Byiii, all other 
class A Bianchi models have fractional Kovalevskaya exponents. 



• With non stiff matter (0 < F < 2) class A Bianchi models have at least 
one real or complex Kovalevskaya exponent, except for Bi with fractional 



values of F and Bn with F fractional in 



IH 



/73 



0.82 



These integrability indications arc conforted by (12, for Byiii using Morales- 
Ruis Theory and by using Painleve's analysis for Bjx- 



4 Exact solutions 



4.1 Bi dynamics : The fundamental state 

In vacuum, the general solution of Bj model could be explicited. The line 
element is 

ds^ = t^Pi dxl + t^P^ dxl + i^Pa ^2,2 _ ^^2 ^^^5-) 

where 

{pi = -nu/{l + u + u'^) 
P2 = ^{l + u)/{l + u + u'^) (16) 
P3 = f7u(l +u)/(l +1t + u2) 

The real is directly proportional to Universe Volume variation wich is constant 
in Bj. Toward singularity [t — s- 0), xi expands when X2 and x^ contracts, all 
with a constant exponential rate. Any couple [m, il] G [1, +oo[xR associated to 
an order of axis defines a Kasner state. 



4.2 Bii dynamics : One Kasner transition 



As noted by (3) , Bn dynamics corresponds to a transition between 2 asymptotic 
Kasner states : 



[u,n] (DAO) 

Pi < P2 < P3 

□ A 
Initial Kasner State : t — s- +oo 



[u - 1, n{l - 2pi)] (AOD) if M > 2 

{u - 1)"^ , n{l - 2pi)j (ADO) if u < 2 

Final Kasner State : t — > 
(17) 



4.3 Conjectures (partially proven) 

As indicated by the rigorous proof by JHI ; sll classes of Bianchi models converge 
generically toward a simple Kasner state or a closure of Kasner states when 
< — > 0. Numerical analysis, more precisely Billiard analogy (see next section), 
let us think that : 
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• -Bvio and -Byiio dynamics correspond to a finite number of Kasner transi- 
tions 

• ^viii and Bix dynamics correspond to an infinite number of Kasner tran- 
sitions 



5 Billiard analogy 

Those works was pionneered by 0], and more recently by |13| or |15j . Intro- 
ducing the super-time i such that di = V^^'^dt, a "mass" m — V^l'^ and an 
"energy" E — [dV / dt)"^ /2V'^/'^ , dynamical equation © in vacuum becomes 



dqi 



dt 



dpi 



dt 
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dE 

dE 
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with 



E = 



2m 



1=1 

normal projector onto (e i, 62) plane 

(18) 



Toward singularity, E —> +00 and m — 0. When we go through time backward, 
equations (|18() are ones of a ball with a decreasing mass m and with an increasing 
energy. The exponential nature of the potential ^ allows a precise description 
of he dynamics, namely, the billiard analogy. When ^ is exponentially negli- 
gible, solutions of lfTH|l are almost straight fines in phase plane (51,92) : these 
are almost Kasner states in BKL representation. When ^ cannot be neglected 
anymore, it is mainly due to one of it exponential term, equation of motion — 
for a linear combination y of (51, (72) — takes a form equivalent to^ 



'^^y J.2 y 



with y (0) = — ^ 
dx 







x=0 



This equation could be solved analytically in 



y (x) = In 



= -2 In 



(19) 



(20) 



It shows clearly a transition between the two Kasner states associated to the 
straight lines which are the asymptotes of this solution (see Fig. |2J|. This 
generalized Bu transition corresponds to a bounce of a ball in an amazing 
billiard. As a matter of fact, when p\+p^ — 0, all ball's energy is concentrated 
in the potential term ^, this situation corresponds to the maximum of the y 
function plotted in Fig. El : this is a bounce. It happens effectivelly against 
the isocontour line ^ = C(q*) where dq*/dt = 0. As we go toward singularity 
{t ^ 0) these lines move outward as one can see in Fig. O The Universe ball 
then moves in a expanding billiard. It is easy to prove that ball's velocity norm 
I d(\/dt I is always greater than the one of cushions which is the energy variation 



The X variable represents the supertime t such that when t ^ x one has § ■ 
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Figure 2: Reflexion on a straight cushion (left) and in a corner (right) 



(see |18) for instance). This property remains true when barotropic matter fills 
Universe provided that F < 2. The amplitude of the bounce depends on k : 
the largest k is, the smallest is the bounce angle between the two asymptotic 
linear regimes. Between the two ideal cases (quasi Kasner state and bounce) the 
parameter k is dynamical, its variations remains negligible during each phase. 

6 Dynamical properties of Bianchi Billiards 

In the billiard formulation, dynamics of Bianchi Universe is fully understood 
considering isocontours of the potential ^ which are presented on Fig. |3 Excluding 
Bix and Byiii a-U isocontours are open curves. Hence after a finite number of 
bounces, the ball representing Universe find a hole and go away to reach the 
asymtotic regime solution of equation H19|) which is a Kasner state. By oppo- 
site, for the two closed curves potential, there is no hole in the billiard, cushion 
allways expands less quickly than the ball, and there is an infinite sequence of 
Bii transitions toward singularity. Although -Byiii cushions seems open along a 
vertical channel, the bounce process makes the ball which enter the channel goes 
back. As sketched on Fig. |21 in the channel fc ~ (72, then as the ball sinks into 
this region, (72 grows and bounces are more and more pinched. By consequence, 
as the cushion is not strictly vertical, the ball finally goes back and leaves the 
channel. These results are surely well known and one remaining question could 
be : Is -Byiii dynamics equivalent in some precise sence to Bix one ? In terms 
of complexity the two dynamics seems equivalent as they have exactly the same 
set of Kovalevskaya exponents. Let us show that they are also equivalent by the 
fact that they present the same kind of transition to chaos, and they have the 
same kind of attractor. 
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Figure 3: Isocontours of potential ^ for all distincts Bianchi Universe in vacuum. 
Arrows indicate the increasing values of ^ 
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6.1 Bix and -Bviii Poincare map analysis 



Detained analysis show that Bianchi Universe volume is a regular variable. It 
is the only variable attached to cushions expansion in the billiard formalism. 
Therefore, the complexity of the dynamics is contained in the cushions's shape. 
In order to study this complexity in the two relevant cases which are Bix and 
-Bviii Universes, we have studied dynamical propertied of temporal sections of 
the hamiltonian E involving in equation H18|) . As V, t, i and r are all related in 
a bijective way, by temporal section we mean sections of the dynamics described 
by equation H18|l with m = est and E = est. For convenience we use always 
a ball with unit mass, we focuss attention under Poincare maps associated to 
growing values of E. Such maps are obtained by considering the intersection 
between an orbit (a solution of the 4D differential system) and a fixed plane, 
namely 11 — {q2 — 0,p2 ~ 0} for maps plotted in Fig. ^ The greater are values 




Figure 4: Poincare maps toward singularity of temporal sections of to = est and 
E — est of Bix (top : TO = 1 and E = —2.8, —2.1, from left to right) and Bviii 
(bottom : m — 1 and E — 5.0, 10.0 from left to right) dynamics. For each map 
a set of 9 homogeneously distributed initial conditions have been choosen. 

of the energy E, the closest from singularity sections are. Analysis of such maps 
is clear using KAM Theory : When periodic or quasiperiodic orbits correspond 
to Poincare sections disposed along curves, chaotic orbits sections fill dense 
regions. We can then conclude that far from singularity orbits are regular in 
both Bix and i?viii- Approaching singularity, a transition to chaos seems to 
appear in the dynamics of both these closed cushions expanding billiards. In 
order to produce a comparative measure of this chaos we propose to produce 
an analysis base on Bix and -Byiii Universes truncated dynamics. 
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Figure 5: Initial conditions 9 and lo for fractals maps represented on Fig. 
6.2 Bix and -Bviii fractal attractors 

As suggered by all previous indicators (Kovalevskaya exponents, Poincare maps), 
we suspect i?viii dynamics to possess the same kind of attractor exhibited ini- 
tially by |17j for i?ix- In order to confirm this suspition we have applied the 
truncated dynamical technique used by in the context of Hamiltonian for- 
malism described above. 

Initial conditions we used are = [O, — V61n 2, 51] and Po = [cos 6*, sin 6*, cj] 
where the constraint H = is fuUfilled provided that 



The upper and lower signs are associated respectively to Bix and Bym Uni- 
verses. The two free parameters uj and 9 fixe respectively the initial variation 
rate of Univers volume (or the initial escape speed of cushions in the billiard 
formalism) and the orientation of the initial velocity vector. These initial condi- 
tions are represented on Fig. El The numerical integration of the dynamics from 
these initial conditions shows an unlimited sequence of Bu transitions between 
Kasner states (see formula (|17ll '). We stop numerical integration when the Kas- 
nerian u parameter is greater than a critical value (for maps presented on Fig. 
El we have used u^xu = 8). When the dynamic is stopped, we affect to the corre- 
sponding point in the 9 — uj plane a color which correspond to the expanding axe 
of stopped Kasner state (Red <-> Ai, Blue ^ A2 and Green ^ ^3). We have 
studied the [0, 7r/2] x [—2, —3] portion of the 9 — u; plane, projected on a regular 
500 X 500 grid. The corresponding maps are presented in Fig. El From this map 
we can compute the 3 Haussdorf dimensions associated to each color set. We 
can resume these three numbers into one by averaging each dimension weighted 
by the surface proportion occupied by the corresponding color. We then obtain 
a mean Haussdorf dimension of each map which is d = 1.6976 ± 8.9 x 10^^ for 
i?viii's map and d = 1.7141 ± 8.5 x 10^'^ for i?ix's map. These two numbers 
are then numerically indistinguable. In this sense we claim that -Bym and Bix 




(21) 
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Figure 6: -Bviii (left) and Six (right) fractal maps in the [0,7r/2] x [—2,-3] 
portion of the 6 — w plane 

universes are dynamically equivalent 
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